Landau-Zener-Stueckelberg effect in a model of interacting tunneling systems 
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The Landau-Zener-Stueckelberg (LZS) effect in a model system of interacting tunneling particles 
is studied numerically and analytically. Each of TV tunneling particles interacts with each of the 
others with the same coupling J. This problem maps onto that of the LZS effect for a large spin 
S = TV/2. The mean-field limit TV — > oo corresponds to the classical limit S — * oo for the effective 
spin. It is shown that the ferromagnetic coupling J > tends to suppress the LZS transitions. For 
TV — > oo there is a critical value of J above which the staying probability P does not go to zero in 
the slow sweep limit, unlike the standard LZS effect. In the same limit for J > LZS transitions 
are boosted and P — for a set of finite values of the sweep rate. Various limiting cases such as 
strong and weak interaction, slow and fast sweep are considered analytically. It is shown that the 
mean-field approach works well for arbitrary TV if the interaction J is weak. 

PACS numbers: 03.65.-w, 75.10.Jm 



I. INTRODUCTION 

The Landau-Zener-Stueckelberg (LZS) problemi*2^ 
of quantum-mechanical transitions at an avoided level 
crossing induced by a linear-in-time energy sweep is well 
known in many areas of physics, in particular, in the 
hysics of atomic and molecular collisions (see, e.g., Ref. 
Recently the LZS effect was observed in crystals of 
single-molecule magnets and it was used to extract their 
tunneling level splittitng 

The specifics of the LZS effect in solid-state systems 
is the macroscopically large number TV ot tunneling 
species such as spins of magnetic molecules in the above 
experiments^^ These tunneling species are usually cou- 
pled with each other (mainly via the dipole-dipole in- 
teraction in the case of magnetic molecules) and the 
coupling energy J can easily exceed the splitting A. In 
this situation tunneling species (that can be described 
by pseudospins S = 1/2) do not tunnel independently, 
and one has to consider the Schrodinger equation for 
the whole system that is described by 2^ variables. 
The latter is a new and tremendous problem since the 
whole Schrodinger equation becomes intractable even for 
a moderate number TV while analytical treatment is dif- 
ficult because the system is far from the equilibrium. 

As a plausible first step towards the solution of the 
LZS problem for a system of interacting particles one 
can consider a simplified model in which each particle 
interacts with all TV — 1 other particles with the same 
strength J, as was done in Ref. For this model the 
Schrodinger equation simplifies so that one has to solve 
only TV + 1 equations instead of the 2^ equations in the 
general case. In the limit TV — > oo the problem simpli- 
fies since the mean-filed approximation (MFA) becomes 
exact. The latter allows one to gain insights into the 
problem by analyzing its mean-field solution numerically 
as well as analytically in different limiting cases. Because 
the molecular field exerted on a tunneling particle from 
the others depends nonlinearly on time in the region of 
the level crossing, the problem cannot be linearized (c.f. 
Ref. 0) and the results differ essentially from the stan- 
dard LZS solution. For instance, the ferromagnetic cou- 



pling, J > 0, makes the energy sweep faster, and thus the 
probability P to stay on the initial bare level increases. 
The nonlinear time dependence of the sweep within the 
MFA has stimulated Ref. |8|, where the different kinds of 
sweep nonlinearities were investigated for a single tunnel- 
ing system. Certainly the MFA for systems of interacting 
tunneling particles is a more complicated problem since 
in this case the form of the sweep is unknown from the 
beginning and it should be found self-consistently. 

It is very interesting to study the difference between 
the exact quantum-mechanical and the mean-field solu- 
tions of the LZS problem for finite TV. The original model 
maps of the model of a large spin S — TV/2 with the uni- 
axial anisotropy D ~ J , transverse field H x ~ A, and the 
sweeped longitudinal field H z (t). That is, the mean- field 
limit TV — > oo of our model corresponds to the classical 
limit S — ► oo for the effective large spin. In Ref. it 
was shown numerically that for the model with TV = 4 
and J > the exact and MFA results are qualitatively 
similar. It is very interesting, however, to investigate the 
problem for larger TV and for J < 0. 

Numerical and analytical study of these problems is 
the objective of this paper, the rest of which is organized 
as follows. In Sec. [H] the Hamiltonian is written down, 
the Schrodinger equation is simplified making use of the 
symmetry of the interaction, the problem is mapped onto 
that of a large spin and its general properties are stud- 
ied. In Sec. II I II the case of well-separated resonances, 
| J| ^> A where the problem can be reduced to that of suc- 
cessive standard LZS transitions is investigated. In Sec. 
IIVI the opposite case J <C A, where the problem can be 
solved perturbatively in J, is considered. Final results are 
worked out in the cases of fast sweep and slow sweep. In 
Sec. [V] general properties of the mean-field solution of the 
LZS problem with interaction are investigated, in partic- 
ular, with the help of the mapping onto the classical-spin 
problem. Numerical solution shows, in particular, that 
for J < complete LZS transition, P = 0, is achieved at 
some values of the sweep rate. Analytical treatment of 
the slow-sweep limit within the MFA is provided in Sec. 
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II. THE HAMILTONIAN 



for k = 0, . . . , N and with the bare energies 



We consider the model of N double-level tunneling sys- 
tems described by pseudospins 1/2 and interacting each 
with each with an equal strength 

1 N J N 

H = -^[W(t)a iz + Aa ix ]- - £ a iz a jz . (1) 



i=l 



Here er^ are Pauli matrices, W(t) is the energy sweep that 
is taken to be linear in time 



W(t) = -£i(i) = vt, 



(2) 



A is the level splitting at resonance in the absence of in- 
teraction (t = and J = 0), whereas J is the interaction 
constant. The case J > corresponds to the ferromag- 
netic (FM) coupling, whereas that of J < corresponds 
to the antiferromagnetic frustrating (AFMF) coupling. 
In the latter case, the ground state of the system with 
W = A = is a highly degenerate state with the minimal 
total spin. If J = 0, the problem simplifies to the well 
known LZS problem for individual tunneling systems. If 
these system at t — > — oo are in the bare ground state 
= | i) before crossing the resonance, the probabil- 
ity to stay in this state after crossing the resonance at 
t — > oo is given by,2ii2iiS 



P(oo) = P = e 



ttA 2 
2hv 



(3) 



(Probabilities without the argument are shortcuts for the 
final-state probabilities at t — oo throughout the paper). 
For the problem with interaction, the wave function of 
the system can be written as the expansion over the 
direct-product states 

mi ,...,mjv=— 1)1 



i = (J) = 11), tfi = (J) =IT>- (4) 



The initial condition for \&(t) is C_i l ... l _i(— oo) = 1 
whereas all other coefficients are zero, i.e., the system 
starts in the state with all spins down. With time the 
state of the system becomes a superposition of all possi- 
ble basis states in Eq. The one-particle probability 
to remain in the initial state —1 for our iV-particle system 
is given by 



P N (t) 



E 



|C-i, 



(5) 



,mjv — — 1)1 



and it starts from Pn(~oo) = 1. 

The solution of the LZS problem for this model is sim- 
plified by the fact that the coefficients C mil ... )Tnjv (t) de- 
pend only on the number k of spins up while they are 
independent on the choice of these spins. Thus one can 
label the states by the index k only that results in the 
Schrodinger equation 



ihC k =E k (t)C k 



kA 



-Q 



k-l 



(N - k)A 



C, 



k+l 



(6) 



E k (t) =f^-kj W{t) + 2Jk(N -k) + const. (7) 

Each fc-state is realized by Nl/ [(N — fc)!fe!] different 
choices of the indices mi, . . . , tun in Eq. Q}. For J < 0, 
the ground state with W = A = corresponds to 
k = N/2 and is highly degenerate. In Eq. © mi is 
fixed to —1, and the remaining indices ma, . . . , mjv can 
also be parametrized by k. In the corresponding number 
of realizations, one should use N — 1 instead of N that 
leads to 



N-l 



N 



(') = £ 



(N-iy. 



fc=0 



(N-l- k)\k\ 



\c k (t)\' 



It is convenient to introduce the coefficients 



Nl 



('k 



(N — k)\k\ 



C k 



(8) 



(9) 



that satisfy the normalization condition 1 = X^fcLo l c fc| 2 
and the Schrodinger equation 

ihc k = E k (t)c k - — Zfc-i,fcCfc_i - yZfc, fc+ ic fc+1 (10) 



with l k)k+ i ee y/(N - k)(k + 1). Then Eq. © simplifies 
to 



fc=0 



k 



P N (t) = ( 1 - Jj ) Pk(t), Pk(t) EE \c k (t)\ 2 . (11) 



Eq. (|10|) maps on the corresponding Schrodinger equation 
for a large pseudospin S — N/2 with the Hamiltonian 



H s = -H z {t)S z - H X S X - DS 2 Z + const, (12) 



where 



H z (t) = W(t), 



H , 



A, 



D = 2 J 



(13) 



(we set gfi B — 1 so that the "magnetic field" is energy 
dimensional). Note that the excitation number k is re- 
lated to the eigenvalue m of S z by k — N/2 + m. Thus 
P(t) maps onto 



(14) 



The large-spin model defined by Eq. I|12l) is well known 
from the physics of the single-molecule magnet Mn 12 . For 
small A the energy levels of the system as function of W 
are nearly straight lines with avoided crossings (see Fig. 
0). Splittings at avoided crossings are or order A(A/J) ra , 
where n are appropriate powers^! 

Our model described by Eqs. or l|12(l . although not 
completely realistic, allows one to trace down the influ- 
ence of the interaction on the transition probabilities and 
to obtain a number of interesting analytical results. We 



3 



k=2 


Ferromagnetic coupling </ ^/ 
\ N=4, A = 0.l /y^ r 


a 






/V\3 








k =Y \k = N=A 




0.6- 



0.4 



0.2 



A= 1 
iV=3 
/ = 10 



B 

o 



5 



-4-3-2-10123 .„ 4 




FIG. 1: a - Exchange-split resonances for a ferromagnetically 
coupled spin cluster for |J| > A. Transitions shown by the 
dashed and dotted arrows arise in higher orders in A/ J and 
are much weaker. 6 - Same for an AFMF coupled spin cluster. 



will see that, in particular, suppression of transitions for 
the ferromagnetic interactions (J > 0) can be easily un- 
derstood. An advantage of this simplified model is the 
possibility to solve the problem numerically up to much 
larger system sizes, compared to models with realistic 
interactions. For the latter, the number of different co- 
efficients in Eq. Q and thus of differential equations to 
solve is 2^, whereas in our case there are only N + 1 
equations. 

For numerical calculations we use Wolfram Mathemat- 
ica 4.0 that employs, in particular, a very accurate differ- 
ential equation solver needed for handling strongly oscil- 
lating solutions over large time intervals. The results of 
solving Eq. I|1U|I are shown in Fig.[21for a cluster of three 
ferromagnetically coupled particles. 



III. WELL-SEPARATED RESONANCES 

The interaction can profoundly influence the LZS ef- 
fect if it is strong enough, | J| > A. The general tendency 
can be already seen from the well-known mean-field ar- 
gument. If one of the tunneling particles flips to another 
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FIG. 2: Time dependence of the one-particle staying proba- 
bility Pjv as function of the energy sweep W(t) for a system 
of N — 3 ferromagnetically coupled particles. At intermedi- 
ate sweep rate, e = 1 only the first-order LZS transition at 
W — 2Jo (the A transition) is seen. For rather slow sweep 
rates (here e = 10), higher-order transitions (here the A 2 /J 
transition at W = 0) are switching on. 

bare state during the resonance crossing, then the total 
field (the external sweep field plus the molecular field) on 
all other particles jumps. For the ferromagnetic coupling, 
J > 0, the jump of the total field is positive, other parti- 
cles are brought past the resonance and lose their chance 
to flip, and thus transition probability is suppressed. For 
the antifcrromagnetic frustrating coupling, the jump of 
the total field is negative and other particles are being 
set back before the resonance and are getting one more 
chance to cross the resonance at some larger field value 
and flip; Thus the transition probability for the system 
should increase. 

The influence of the interaction can be readily illus- 
trated within a rigorous many-body quantum language 
for our model if one considers the energy levels of the 
system shown in Fig.^ One can see that instead of a sin- 
gle resonance at W = for individual or noninteracting 
particles there is an interaction-split resonance that con- 
sists of many avoided line crossings. These avoided line 
crossings are well separated from each other on the plot 
if | J | ^> A. The question whether these well-separated 
crossings can be considered as a succession of indepen- 
dent LZS transitions (i.e., whether they are dynamically 
well separated) depends on the sweep rate. Whereas for 
the slow sweep, e > 1 [see Eq. (0] the condition \J\^>A 
is sufficient, for the fast sweep, e < 1, a more stringent 
condition is required that follows from Eq. (16) of Ref. 
8. The resulting combined condition for the dynamically 
well-separated resonances thus would be 

|J|>> {^A/e, ell. ^ 

This is, however, only an apriori estimation considering 
two resonances. For N ^> 1 there are many resonances, 
and the deviations from the single-resonance picture can 
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accumulate with the increase of N. One can do an aposte- 
riori estimation in the slow-sweep limit where the result 
for the well-separated resonances P — > 1 — 1/ N for e — ► oo 
following from Eq. (|2U|) as well as the mean-field result 
for N ^S> 1 and non-separated resonances [the first line 
of Eq. are available, in both cases 1 — P <C 1. The 
crossover between these results should take place at |J| ~ 
AN 1 / 2 , thus the first line of Eq. lj*Po)l should be replaced 
by 



I J| » AN 1 / 2 , 



e > 1. 



(16) 



The second line of Eq. i|15|) could be also modified by N 
but it is difficult to derive an exact condition. 

A perturbative analysis of the stationary Schrodinger 
equation for our model shows that tunnel level splitting 
of the bare levels with different values of k has the form 11 



SE, 



fei .k,2 



J 



1*1— &a|— 1 



(17) 



Thus for | J | ^> A only the first-order or direct transitions 
with \k\ — fol = 1 should be taken into account whereas 
the higher-order transitions are weak (see Fig. . 

Now it becomes clear that for the ferromagnetic cou- 
pling and well-separated resonances only the strong tran- 
sition between the initial level k = and the level k = 1 
happens. That is, in another language, only one particle 
of N has a chance to tunnel, and the tunneling proba- 
bility for the whole system is strongly reduced. For a 
quantitative analysis one can neglect all c k with k > 1 in 
Eq. I|1(J|) that yields the system of equations 



ihcn 



E c 



A 



ihci = E\C\ — —VNcq 



(18) 



that maps on the standard LZS problem with A => An = 
Ay/N and with the resonance at W = 2(N — I) J instead 
of W = 0. The final-state probabilities for Eq. lfT%|) are 
according to Eq. © 



P, = P N = e~ N °, 



Pl = l-P 



N 



(19) 



Then with the help of Eq. (|1 1|1 the one-particle staying 
probability for J > can be cast into the form 



P. 



N 



1 



N 



[1 - P N ) 



(20) 



and it varies between 1 in the fast-sweep limit and 1 — 1/ N 
in the slow-sweep limit. Thus suppression of transitions 
by the ferromagnetic coupling is extremely strong for a 
large number of particles N. Expansion of the transition 
probability of Eq. (|20|l in the fast-sweep limit reads 



1 



Pn = e-^e 2 + 0(e 3 ), 



(e < 1/N). (21) 



Note that the first term of this expansion is insensitive 
to the interaction (in this context to the number of in- 
teracting particles) and is the same as for the standard 



LZS problem, 1 - P = 1 - e 



72! + e 3 /3! 



n= io,y>o 



Separated resonances 




FIG. 3: One-particle staying probability Pn vs the sweep- 
rate parameter e for systems with FM and AFMF interaction 
in the case of well separated resonances. 



The AFMF coupling J < favors transitions, as can 
be seen from Fig. The original resonance is splitted 
by the interaction into N resonances filling equidistantly 
the range -2(N — 1) J < W < 2(N - 1) J. One can see 
that in the limit of slow sweep the system remains on the 
lowest adiabatic energy level thus Pn —> 0, in contrast 
to the ferromagnetic case. For well-separated resonances 
the problem described by Eq. llUfl splits into indepen- 
dent standard LZS problems for the resonances between 
the levels k and k + 1 that are described by an effective 
splitting Ajv,fc = AZfe^+i [cf. Eq. (jT%|l ]. It is convenient 
to designate the probability to stay in state k after cross- 
ing with state k — 1 but before crossing with state k + 1 
(see Fig.njj) aspfc(mid). Then solutions of the LZS prob- 
lems for individual crossings along with conditions of the 
probability conservation can be written as 

Pk (oo) = P( jv - fe )( fc+1 )p fe (mid) 
Pfc(mid) = p fc _i(mid) -p fc _i(oo) 
Pk -i(oc) = P^-^V-iCmid), 

etc. These can be combined into the recurrence relation 



(22) 



Pk 



(oo) = P^-k)(k+i) (p-W-k+Dk_ 1 j pk _ 1 ( oo) 
= P w - 2fc (l-p( w - fc+1 ) fe )pn(co) (23) 

that has to be iterated with the initial condition po(°°) — 
P N . In the slow-sweep limit e 3> 1 /N, one has P N <C 1 , 
thus one can drop the factor (l — p( N - k + x ) k ^ ^ l in the 
second line of Eq. ll'MI) that after iteration results in p k = 

p{N - k) (k+i)^ In this cas6i in Eq (jnji p ^ pn-1 = p n 

are dominant for large N whereas all other summands 
are much smaller. This yields 



Pn = 



1 

1+ N 



P N « 1, 



(e > 1/JV) 



(24) 



for the antiferromagnetic coupling. In the fast-sweep 
limit the solution of Eq. (|23|l can be expanded in pow- 
ers of e, after which Eq. I|ll|) yields for the transition 
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FIG. 4: Numerically computed staying probability Pjv vs 
J/ A for N = 3 and different sweep rate parameters e. The 
horizontal dashed lines on the left side of the plot are asymp- 
totes corresponding to well-separated resonances. 



probability 



37V 



l-PAr = e+^-2je 2 + 0(e 3 ), (e«l/A0- (25) 

Again, the first term of this expansion is the same as 
in the standard LZS problem. For intermediate sweep 
rates one can iterate Eq. ill'-il) numerically and substitute 
the solution for pk into Eq. The resulting curves 

Pjv vs e are shown in Fig. [3] along with those for the 
ferromagnetic coupling, Eq. (|20|l . 

To conclude this section, we show the numerically com- 
puted dependences of Pjv on J/ A for A = 3 and different 
e in Fig. 0] The case of well-separated resonances that 
was considered in this section corresponds to the plateaus 
on the left and on the right sides of the plot. 



IV. FAST SWEEP AND WEAK COUPLING 

In the preceding section we have considered the strong- 
coupling limit | J| ^> max (A7V 1 / 2 , A/y^e) in which indi- 
vidual resonances are well separated from each other and 
one thus deals with successive standard LZS transitions. 
The opposite limiting case is that of the weak coupling J. 
This case can be solved by a perturbative expansion in J 
with the standard LZS solution as the zeroth-order ap- 
proximation. The situation further simplifies in the fast- 
sweep limit, where the transition probability is small and 
the coefficients Ck in Eq. ljTU|> decrease with k as powers 
of e <C 1. In particular, cq contains terms of orders e°, e 1 , 
and e 2 , the coefficient c\ contains e 1 / 2 and e 3 / 2 , the coef- 
ficient C2 starts with e, etc. Knowing these contributions 
into Co, ci, and C2 is sufficient to set up the expansion of 
Pjv up to the first nontrivial order e 2 . A straighforward 
but cumbersome calculation yields 



1-P 



N 



4Jo 
\J1~kYlv 



0(£ 3 ) 



4J 
AT' 



.5/2 



where we have defined 



Jo = (A - 1) J. 



(26) 



(27) 



One can see, again that ferromagnetic interaction sup- 
presses transitions whereas the AFMF interaction facil- 
itates transitions, and that the effect of interaction is 
increased by the number particles in the system. In 
fact, however, in our model the interaction should scale, 
on physical grounds, with the system's size, i.e., Jo = 
(N — 1) J should be size independent. Note that in the 
fast-sweep limit one could do the expansion in e for ar- 
bitrary J. This leads, however, to cumbersome multiple 
integrals. The strong-coupling limit for the fast sweep 
e«1 has been considered above, and the corresponding 
counterparts of Eq. I|26|l for ferro-and antiferromagnetic 
coupling are Eqs. (|21() and l|25|) . 

It should be noted that although Eq. 1|26[) is valid for 
e<l, the method of its derivation above is only justified 
for £ <C 1/A, i.e. for much faster sweep rates, if N ^> 1. 
Indeed, in the course of the derivation the terms of orders 



(eN) 



1/2 



eN, etc., appear that are required to be small. 



Only at the very end the leading N contributions cancel 
each other that leads to Eq. 12tj|) that fortunately has a 
larger applicability range. In fact, for large systems with 
weak interaction the most populated final states are 



= A(l - P) = N£, 



(28) 



so that keeping only the states with k = 0, 1, 2 above was 
not justified for Ae > 1. Eq. I|28|l can be easily derived for 
an assembly of noninteracting tunneling species. To this 
end, we use the coefficients Cf. of Eq. © and represent 
them in the factorized form 



C k 



N-kk 



1 ' 



(29) 



where a_i and a\ arc the coefficients of the wave function 
of the one-particle problem, ip = a-\ip_ 1 + a\^ 1 . This 
yields 



|C fc l 2 = 



-i| 2(JV - fc) |ai| 2fc =P"- fe ' 



c (l -P) k 

that can be plugged into Eq. JSJ to give 
(A-l)! 



Pv = 



JV-l 

^ (A- 



1 - k)\k\ 



pN-k (1 



P) 



(30) 



(31) 



The summand in this formula has its maximum at k = 
fcmax given by Eq. (|2*S|l and it becomes a narrow Gaussian 
packet around fc max for N 1. 

A more rigorous method of handling the fast-sweep and 
weak coupling limits uses the one-particle density matrix 
that is defined by 



= P 



P-i,-i 
with Pjv of Eq 

Pi,-i = 



N- 



Pl.l = 1 

| and 



p-i.-i 



l-P 



(32) 



m2,...,m N 



6 



N-l 

E 

fc=0 

pi,-i 



(N-iy. 

(N-l- k)\k\ 



(33) 



The density-matrix equation (DME) can be obtained by 
time differentiating of p rnn and employing Eq. @ that 
yields 



h P-i,-i 



iA 



A Im p 1 



ft'Pi-i = ~iW(t)p 1 _! + — ( Pl l - p_ x _i) 



2iJ 



(N-l)\ 



N-l 

E '(jv- 1 - fc)!fc! 



(N — 2k - 1) C fe C fe * + ^34) 



In the absence of interaction, J = 0, the last term in 
p_l i disappears and one obtains a DME for one isolated 
particle that is equivalent to the one-particle Schrodinger 
equation. Note that solving this one-particle equation in 
the fast-sweep limit requires e <C 1, in contrast to Eq. © 
with J = that requires e <C 1/N. 

For J ^ Eqs. (|34[1 do not form a closed system of 
equations. In this case Eqs. (|34[1 can only be useful if 
an approximation for the interaction term be found. In 
particular, one can consider the interaction term pertur- 
batively in J and use 



p = 

rrnn rmn 



S Pr, 



(35) 



where pmn is the density matrix without interaction that 
satisfies 



P-1,-1 



(36) 



and <5/o m „ is the correction term. The latter satisfies 
equations similar to Eqs. (|34|l in which, however, CkC^ +1 
is replaced by its noninteracting-particle expression fol- 
lowing from Eq. 



a N 7 k a* (a* 



, N-k-l 



\k+l 



I |2(JV — fc — 1) I \2k * 

= \a-i\ y >\ai\ a_ia x 



- I o {0) 



( a l) 

* 
1 

(l-p%,-i)"p%- (37) 



With the use of the latter, the sum over k in Eqs. 134|) 
can be performed to yield 



h5p_i_i = AlmSp 1 _ ± 

Mp 1 _ x = — iW(t)Sp_i x — iA8p_ x —1 

-2ijJl-2 P ^_ 1 )pfl v 



(38) 



Solving Eqs. I|3t)l) and (|38|) perturbatively for e <C 1 re- 
sults in Eq. igSJ. 

In accord with the remark at the beginning of this sec- 
tion, J-perturbative Eqs. (|36l) and H38|) can be solved in 
terms of the hypergeometric functions for any e (For the 
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FIG. 5: Numerically computed staying probability Pjv vs 
p = 2 Jo/ A for e = 1 and different system sizes N, including 
the mean-field result for N — > oo. 



standard LZS problem, J = 0, this was done by Zener, 2 ) 
It can be shown that in the slow-sweep limit this solution 
simplifies to 



dP 
dp 



p=0 



2 J 

cee- £ , P=~^> ( 39 ) 



where c is a number and the parameter p should not be 
confused with the density matrix. 



V. THE MEAN FIELD APPROXIMATION 

If each of tunneling particles interacts with all the 
other particles with a coupling of the same sign, as is 
the case in our model for N 3> 1, then the mean- field ap- 
proximation (MFA) should work well. The MFA employs 
a one-particle density-matrix equation or a one-particle 
Schrodinger equation in which the interaction enters via 
the molecular field. For the model of Eq. (JTJ one has 
W{t) =*> W cS (t) where 

W cS (t) = W(t) + 2J (l-2p_ 1 _ 1 ) (40) 

in the DME of Eq. HJ, i.e., 

h P-i-i = Almp 1 _ 1 



h'Pi.-i = -iWafftfJp!,.! + iA I - - p_ lt _ x j .(41) 
This is equivalent to the nonlinear Schrodinger equation 
ihC-i = iWe ff (t)C_i - jCi 



1 



-C-l 



(42) 



with W e g(t) = W{t)+2J {l-2\C- 1 \ 2 ). It is interesting to 
note that if one solves Eq. I|41|) perturbatively in J using 
Eq. (1231 one obtains Eqs. (O and ||3BJ. This implies 
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that the MFA works well in the weak-coupling limit in 
our model (see Fig. [SJ) . 

In the case N 3> 1 one can assume that the state of the 
system is described by a narrow packet in the fc-space, 
as is indeed the case for noninteracting particles, see Eq. 
(|3l"1l . Then in Eq. one can replace (N - 2k - 1) => 
(N — 2fc max — 1) with k max of Eq. I|28|) . after which the 
sum over k converts to p_ 1 1 according to the definition in 
Eq. (j2SJ). This leads to Eq. ||ITJ. One should, however, 
realize that this is no more than a heuristic derivation 
since it was not proved that the state of a system with 
interaction is described by a narrow packet for N 3> 
1. On the other hand, the narrow-packet assumption is 
quite plausible since the limit N 3> 1 corresponds to the 
quasiclassical limit S ^> 1 of the equivalent model of Eq. 
(|12[1 and the states of quasiclassical systems should be 
narrow packets. It is remarkable that the density-matrix 
equation within the MFA, Eq. (|41|l . can be interpreted as 
a purely classical equation of motion for a spin vector, the 
Landau-Lifschitz equation (LLE). Indeed, rewriting Eq. 
(|33[l in terms of c k with the help of Eq. © and adding 
Eq. i|14|) one obtains the mapping relations 



1.0 



Pi,-i 



1 N 

i = -Yi 

1 N f-i 



k,k+lC k+1 C k 



k=0 



2S 



1 - 



(43) 



Now one can see that the DME in the MFA, Eq. |@TJ is 
equivalent to the LLE for the classical spin vector com- 
ponents 



(44) 



that reads 



HS Z = — ASy 

hs+ = -i [W{t) + 4J(S* - l/2)sj s+ + iAs z . (45) 
The vector form of this LLE is 



s = 7 [s x H cf f] , 



H 



off 



ds 



(46) 



where 7 = 1/h and the effective classical energy is given 
by 



Hcs = ~H z (t)s z - H x s x - D c \s 2 z 



(47) 



with 



H z (t) = W(t), H x = A, 

D cl = D(S - 1/2) = J(N - 1) = J (48) 

[cf. Eqs. H12JI and 1)13(1] . The one-particle staying proba- 
bility Pjy is given by the formula 



(49) 
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that is analogous to Eq. i|14|) . 
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FIG. 6: The mean-field solution for the one-particle staying 
probability P N for the AFMF coupling, Jo = (N — 1) J < 0. 



Let us now analyze the mean-field solution of the LZS 
problem based on Eqs. I|41|) or H46fl and compare it with 
the exact solution of Eq. lfTU|l . We start with the AFMF 
coupling, for convenience. Fig. shows the dependence 
of the one-particle staying probability Pn on the sweep- 
rate parameter e for different coupling strengths Jo . The 
curves on the plot are qualitatively similar to those in the 
case of well-separated resonances, Fig. [3] and they show 
that AFMF interaction boosts transitions, especially for 
slow sweep, e ^> 1. The difference is that in the MFA only 
the combined parameter Jo = (N — 1) J enters, whereas 
for well-separated resonances the results depend on N 
only and not on J. Another difference is that for strong 
enough interaction the mean-field solution for Pn oscil- 
lates and vanishes at some values of e. Oscillations of this 
kind have been found in the solution of the LZS problem 
with nonlinear sweep&iS As was explained in Ref. 0, os- 
cillations arise when the sweep slows down in the vicin- 
ity of the resonance so that (in the extreme case) the 
particle oscillates between the two quantum states and 
the resulting staying probability depends on the phase 
of this oscillation at the end of the stay near the reso- 
nance. This is also the case for the model with AFMF 
interaction in the MFA. Tunneling of particles creates a 
molecular field that changes in the direction opposite to 
the sweep so that the resulting effective sweep W c s(t) 
of Eq. I|40|l can even become nonmonotonic (in this case 
the analytical method of Ref. 7 based on the lineariza- 
tion of the problem near the resonance and introducing 
an effective sweep rate breaks down) . The MFA solution 
becomes time symmetric in the case of complete conver- 
sion, Pn(oo) = that for Jo /A = —10 is realized for 
e = 0.1137983, 0.13207, 0.1741, 0.2233, etc. The cor- 
responding time dependences of P/v(t) and W c s(t) are 
shown in Fig. Q. Although they resemble the complete 
conversion solutions found in Ref. ll3l there are no appar- 
ent analytical solutions for this model. The difficulty of 
the present mean- field model is that W c g (t) is not known 
from the beginning but is a solution of a self-consistent 
nonlinear problem. 
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W(t)/A 

FIG. 7: a - Time dependence of the staying probability Pjv 
for the AFMF coupling J = (N - 1)J = -10 for the first 
three complete conversion cases, Pn (oo) = (see Fig.[§J. b - 
Same for the effective sweep W e s of Eq. 1401 . 



Comparison of the exact and the MFA solutions for 
the one-particle staying probability Pn for the antifcr- 
romagnetic coupling Jo = —10 is shown in Fig. [H] One 
can see that the exact quantum solution converges to 
the classical mean- field solution for N — > oo, although in 
the vicinity of the first complete-conversion point conver- 
gence is extremely slow. 

Let us now turn to the ferromagnetic interaction. The 
mean-field solution for different interaction strengths is 
shown in Fig. [5] Again, qualitatively it is similar to the 
solution in the limit of well-separated resonances, Eq. 
(|20[l that is shown in Fig. [3J The difference is the same 
as for the AFMF coupling: The MFA solution depends 
on J = (N — 1) J whereas Eq. i|2U)l depends on N only. 
Another difference is that in the slow-sweep limit the 
mean- field curve tends to a constant (for J /A > 1/2) 
whereas all quantum curves tend to zero. The latter are 
in fact combinations of several LZS exponentials corre- 
sponding to transitions of different orders shown in Fig. 

by dotted arrows (see also Fig. [2J . One can see, in 
particular, that the curve N = 2 consist of two different 
exponentials, the curve N = 3 consist of three different 
exponentials, etc. With increasing of N, however, the 
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FIG. 8: Comparison of the exact quantum and the MFA 
solutions for the one-particle staying probability Pm for the 
AFMF coupling showing partially slow convergence to the 
classical mean- field limit for N — *■ oo. 
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FIG. 9: The mean-field solution for the one-particle staying 
probability Pm for the ferromagnetic coupling, Jo = (N — 
I) J > 0. The dashed line is the large-e asymptote of Eq. 



coupling J = Jq/(N — 1) decreases and transitions be- 
come not well separated. In this case dependence Pn of 
e becomes a long-tale curve. 



VI. SLOW SWEEP IN THE MFA 

Here we analytically consider the slow-sweep limit 
within the mean-field approximation, to get more insights 
into the mechanism that leads to suppressing transitions 
in the case of FM interactions (see Fig. as well as 
into that of facilitating transitions and probabillity os- 
cillations in the case of the antifcrromagnetic frustrating 
interactions (see Fig. |SJ • 
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FIG. 11: Comparison of the exact quantum and the MFA 
solutions for the one-particle staying probability Pm for the 
ferromagnetic coupling. 



FIG. 10: Comparison of the exact quantum and the MFA 
solutions for the one-particle staying probability Pm for the 
ferromagnetic coupling. 



A. Adiabatic case: Basic equations 

For the one-particle LZS problem, a convenient method 
of treating the slow-sweep limit is that using the adi- 
abatic basis£ Its advantage is that the probability to 
stay on the lowest adiabatic level is at any time close 
to 1, so that one can make an appropriate approxima- 
tion after which the solution is obtained as a quadrature, 
including the case of nonlinear sweep. Direct extention 
of this method for the present model within the MFA 
is cumbersome, however, owing to the complicated self- 
consistent nature of the adiabatic energy levels for a non- 
linear Schrodinger equation or nonlinear density-matrix 
equation, Eq. (|34[) . Fortunately, an alternative descrip- 
tion based on the LLE, Eq. (|46|l . supports a physically 
transparent extension of the method. The adiabatic basis 
for a Schrodinger equation corresponds to the adiabatic 
coordinate system for the classical problem of Eq. 1)46(1. 
The z' axis of this coordinate system is oriented in the 
direction minimizing the classical energy 7i c ff of Eq. I|47|) 
at any fixed time t. It belongs to the x-z plane and it 
makes the angle 6(t) with the z axis that satisfies the 
equation 

h z (t) sin 8 + v sin 8 cos 8 — h x cos 8 = (50) 

following from dH cS /d8 = 0. Here v = Ai/| All = J/\J\ 
and 



K. = 



21 Ad 



A 

2\M 



mo ^ = ™ TO 



2 A 



2 Jo 



Eq. (|50|l is well known in the theory of magnetism. In 
the AFMF case v < 0, its solution is unique. For the fer- 
romagnetic interaction v > 0, Eq. I|50l) has two solutions 
corresponding to the minimum and to the maximum of 
W e ff in the range of reduced fields h x and h z outside the 
Stoner-Wohlfarth astroicU^ (see Fig. 
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FIG. 12: Effective sweep in the adiabatic case for the ferro- 
magnetic coupling of different strengths. 



For h x and h z inside the astroid, Eq. (|50l) has four so- 
lutions corresponding to the stable and metastable min- 
ima as well as to the saddle point and to the maximum 
of 7i e ff • The adiabatic solution corresponds to the mini- 
mum of the energy that changes as h z (t) is sweeped from 
— oo to oo. For h x > 1 the system follows this adiabatic 
solution in the whole range of h z . In this case depen- 
dence of the effective sweep W c g of Eq. (|4"U|) or, here, 
W e e — H z + 2D C \ cos 6, on the energy sweep W is shown 
in Fig. E| One can see that for Jo/ A = 0.5 this de- 
pendence becomes nonanalytic. This corresponds to the 
horizontal line in Fig. 1111 that touches the upper corner 
of the astroid. In contrast, for h x < 1 the adiabatic so- 
lution exists only until the crossing of the right branch 
of the astroid. After that the spin rotates away from the 
disappeared metastable state and its behavior becomes 
nonadiabatic. The latter case will be considered later 
while at first we concentrate on the adiabatic situation. 

The Landau-Lifshitz equation, Eq. 140|) can be rewrit- 
ten in the rotating adiabatic coordinate system as follows 



s' = [s'x ( 7 H^ a + n)] , fl=8e 



(53) 



h 2 J z + h 2 J 3 = 1. 



(52) where 8 is time derivative of the appropriate solution of 
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Eq. 1(50)1 . The reduced effective field 

H e ff 



h e ff 



(h z + vs z ) e z + h x e x , 



2\D cl \ 

in the adiabatic frame is calculated as follows 



(54) 



h^, ff = (h x cos 9 - h e s, z sin 9) e x > 
+ (h x sin 9 + h e g,z cos 9) e Z ' 
h e s,z = h z + v {—s x i sin# + a„> cos9) . (55) 

With the help of Eq. (|50|l components of h' cS simplify to 

h e s,x' = (1 — s z > + s x t tan 9) v sin 6 cos 6 
h c e,z' — h x / sa\9 — s x ivsm9 cos9 

+ (1 - s z ,) {h z - h x cot 9) cos 9. (56) 

Note that if the spin is directed along the z' axis, i.e., 
s x i = s y = and s z > = 1, one has h e g tX r — 0, that 
is consistent with the choice of the adiabatic coordinate 
system. Introducing the dimensionless sweep variables 



W{t) 



vt 
A 



h a 



A 2 

757 < 57 > 



[e should not be confused with e = (w/2)£ of Eq. one 
can rewrite Eq. I(53|) in the form 



ds' 
da 



s' I i 



(58) 



In the slow-sweep limit £> 1 the solution of Eq. I|58|) 
is close to the adiabatic solution s x > — s y — and s z ' = 1. 
Hence one can linearize this equation near the adiabatic 
solution: 



h e «, x ' 
that results in 



1) ft-eff,; 

Sr'V sin 2 



du 

ds y 

du 



sin ( 



h x j sin( 



df9 



(59) 



smf 



where 



C = l-psin 3 6», p = 2J /A. 



(60) 



(61) 



The factor £ in the second of Eqs. I|60|l makes spin pre- 
cession elliptic for any nonzero coupling Jo, especially in 
the extreme case p — 1. This is an important difference 
from the model of one tunneling particle with a nonlinear 
sweep ^ where precession remains always circular. It is 
convenient to introduce 



1 



1/4, 



+ K 



-1/4, 



(l(u) 



and rewrite Eqs. 1(60(1 in the form 

dc+ 1 d9 



du 



-isfl(u)c + 



2 du 



1/4 



(62) 



(63) 



This differential equation can be easily solved with the 
initial condition c + (— oo) = 0, and the staying probabil- 
ity Pjv can be found from Eq. (|49|l . Keeping in mind 
that at t — oo both coordinate systems coincide whereby 
9 = and a = 1, one can rewrite Pjv = P/v(oo) within 
our linearized theory as 



P 



1 



JV 



s 2 ,(oo) +s 2 (oo)l = |£f(oo) 



4 L-x- v--/ ■ 
The final general expression for Pn reads 



P 



- / du^-( 1/4 exr> \ie$(u) 
2j-eo du^ 



where 



du' Q(u ). 



(64) 



(65) 



(66) 



In the absence of interaction J = 0, our solution recov- 
ers that for the standard LZS problem in the slow-sweep 
limit up to the prefactor in front of a small exponential. 
Indeed, in this case the solution of Eq. (|50|l is 



costf = 



y/hj + hi 

in Eq. (|61|l one has p = and £ = 1, hence 

d9 _ 1 
du 1 + u 2 



Q(u) = VI + u 2 , 



2 (u)' 



(67) 



(68) 



One can see that in this case Eq. 1)63(1 coincides with 
second of Eqs. (24) of Ref. and Eq. (|65|l coincides with 
Eq. (26) of Ref. |I with w'(u) = 1 and c_ => 1. The 
latter is exactly the approximation proposed in Ref. 8 
that for the standard LZS problem reproduces the well- 
known result of Eq. J3J for e 3> 1, however with a wrong 
prefactor: P w (7r/3) 2 e~ e . In Ref. we have shown how 
to correct the prefactor by accurately calculating c_ . 



B. Adiabatic case: The results 

In the general case J ^ Eq. 1)50)1 does not yield an 
analytical solution for 9(u). Fortunately, instead of u one 
can use x = cos 9 as the integration variable on the in- 
terval — 1 < x < 1 and with the help of Eq. 1)50)1 express 
u = h z /h x as a function of x. Even better is then to 
parametrize x as x — w/yl + w 2 , analogously to Eq. 
1)67( 1 . This brings Eq. 1)65)1 into the explicit form 



where 



dw 



1 + w 2 



( 1+w 2)3/2 



1/4 



$(wj) 



dwyX 



(1 + w 2 ) 3 / 2 



Ji<S>(w) 



.3/2 



(69) 



(70) 
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In the absence of interaction, p = 0, one has w = u and 
the formulas above describe the standard LZS effect. For 
e>1, the value of Pn is exponentially small and defined 
by the singularities of the integrand in the complex plane 
of w 



Pn = Pme 



where 



4 

Hp) = - 

TT 



dw\fl 



-elm F(p) 



1 - 



(1 



W 2 )3/2 



3/2 



(71) 



(72) 



For p = one obtains the LZS value F(0) = i. For the 
FM coupling the relevant singularity is 



= iy/l -p 2 / 3 , 0<p<l, 



(73) 



that corresponds to vanishing of H(u) of Eq. (|62|l due 
to that of the ellipticity factor £. For p > 1, the motion 
of the spin becomes nonadiabatic (see Fig. Ill|) , and this 
method does not apply any longer. For < p < 1 one 
has KeF(p) = and the limiting forms of ImF(p) are 
given by 



ImF(p) 



2 v / 2 



In- 



:-52 



1 ) p, P < 1 

l-p< 1. 



(74) 



For the AFMF coupling p < 0, there is a pair of relevant 
singularities that also correspond to £ = and are given 
by 



w c± = *(l + |p| 2 / 3 + H 4 / 3 ) 1/4 exp(^ c ) 
1 V3|p| 2 / 3 

= 2 arctan 2TH^- 



(75) 



Note that Im(w c ±) > 1 for p < and thus these singular- 
ities are further from the real axis than the LZS singular- 
ity at w = i. The latter, however, makes no contribution 
since, as can be easily checked, Im F(p) = oo for w c = i 
and p < 0. Analytical calculation of the limiting forms 
of F(p) is more cumbersome for the antiferromagnetic 
coupling. For \p\ -C 1 one obtains 



ReF(p)^2\p\, ImF(p)^l 



, 32 

Vi 1 



(76) 

Note that ImF(p) is given by the same formula for \p\ <C 
1 and both signs of p. In the limit \p\ ^> 1 the result has 
the form 



F(p) 



2 r (V 4 )A„|3/2 



3v^Fr(3/4) 
V^(V3 + i 



,|S/6 



:IpI 



2/3 



r(5/3)r(ii/6)' 

l.l|p| 3/2 - 0.3 |p| 5/6 + (3.6 + 2.1*) |p| 2/3 (77) 
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FIG. 13: Function F(p) of Eqs. (|7T) and l|72|l . 



The real and imaginary parts of F(p) numerically cal- 
culated from Eq. I|72|) . as well as the analytical limiting 
forms, are shown in Fig. 1131 

The prefactor Pnq in Eq. I|71|l is determined for e ^ 1 
by a close vicinity of the singularities w c . For the ferro- 
magnetic coupling the result is 



P, 



NO 



15£pV 1 - P 2/3 



< p < 1. 



(78) 



For the AFMF coupling p < interference of the two 
contributions from u/ c+ and w c - of Eq. (|75|l results in 
the oscillating prefactor 



NO 



I5e\p\\i 



ReF(p) 



p < 0. (79) 



Oscillating prefactors of such a kind take place for one 
tunneling particle if the sweep is nonlinear and decel- 
erating in the vicinity of the resonance^ Here effective 
nonlinearity of this type occurs because of the negative 
coupling between tunneling particles, even for the linear 
sweep. 

Both Eqs. J7S§ and (J7SJ break down in the linit p — > 
since different singularities come close to each other. One 
could work out the crossover from Eqs. (|78() and (|79() to 
the value Pjvo = (tt/3) 2 (see Ref. @) at p = that takes 
place in a narrow region \p\e ~ 1. This makes no sence, 
however, since the result Pno — 7r/3 differs from the 
exact LZS prefactor P — 1 and it has to be improved by 
taking into account nonlinear terms dropped during the 
derivation of Eq. I|tj3|) ■ The corresponding procedure for 
p = is described in Ref. |£| In the present case p ^ 
this would be too involved and we don't try to do it. It is 
clear that the logarithmic singularity of the exponent of 
Eq. H7tjfl should ve compensated for by the singularity of 
the prefactor so that the staying probability Pn behaves 
linearly in p near p — [see Eq. I|39|)]. 

Let us now consider the case 1 — p <C 1 in more detail 
since the prefactor Pjvo 01 Eq. l(78|) diverges at p — * 1. 
Here one has ImF(p) <C 1, and the exponential decrease 
of Pn is very slow. In this region according to Eq. I|73|) 
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the integral in Eq. I|69|l is dominated by small w and it 
can be simplified to 



Pn = U 3/2 



N 

m 

where 

Je1-/)«1 



dt(l + t 2 ) 1/4 e io$ W 
dt'(l + t' 2 f\ 



(80) 



t = 



S 2 ~e. (81) 



It is convenient to compute the integral in Eq. I|8UI) by 
shifting the integration contour by i to suppress oscil- 
lations of the integrand, i.e., to parametrize t = i + z, 
-co < z < oo. With $(z) = z37r/8 this yields 



Pn = P^e-" 11 ^ ^P^voexp 



3tt 



(82) 



where the exponent coincides with the previously ob- 
tained [V3/(2V2)]5 2 e [see second line of Eq. lf73)l]. The 
prefactor P/vo can be calculated analytically for a ^ 1 
and a <C 1. For a» 1 we need the small- z expansions 



6®(z) = $(i + z) - = - (-1 + i) z 5/2 
[l + (z + z) 2 ] 1 / 4 = (-l) 1 /8(2z) 1 /4. 



(83) 



Values of these functions for z < are obtained from 
(-1)5/2 = _■ and (_i)V4 = (l - i)/y/2. After that cal- 
culation in Eq. (|80|l yields 

2nS 3/2 



Pi 



NO 



a > 1 



(84) 



15 a 

that is a limiting form of Eq. I|78|l . In the opposite limit 
a <C 1, the integral in Eq. I|8U|I is dominated by large t, 

so that one can use $(t) = i 4 /4 and (l +i 2 ) 1/4 = i 1/2 - 
This yields 

, , 3tt\ T 2 (3/8) 5 3/2 0.747837 
-'- ] : — 37^' ^ L 



Patq S cos 2 — 



16/ 6^2 a 3 / 4 

(85) 

It is convenient to represent Pno(0'i 8) i n the whole range 
of a with the help of the crossover function /(a) according 
to 



PNo{a,S) 



2tt5 



3/2 



15 a 



-/(«) 



1. 



a > 1 



/(a) ~ i L09294a 1 / 4 , a « 1 



(86) 



Numerically computed function /(a) is shown in Fig. 



C. Strong ferromagnetic interactions 

For stronger ferromagnetic interactions, p > 1, the adi- 
abatic approximation breaks down since the spin approx- 
imately follows the initial energy minimum that becomes 




FIG. 14: Scaling function f(a) of Eq. 



FM coupling, MFA 
slow sweep limit e -> °° 




0.0 0.2 0.4 0.6 



FIG. 15: Pn in the slow-sweep limit e — > oo vs /i^ = 1/p 
A/(2J ). 



unstable at some /i z (t) > (see Fig. 1110 and than it 
performes a large motion that does not approach the 
new stable energy minimum and thus cannot be lin- 
earized around it. While the problem becomes much 
more complicated in this case, one can still find ana- 
lytically lim Pn that is nonzero. To this end, one can 

5 — >OC 

represent the LLE for the spin, Eq. (|46|) . in the Hamil- 
tonian form in terms of the canonical angle variables 
{cos 9, ip} O- {p, q} (in the arbitrary frame) 



• cos 9 



&H d_ 

dt dip ' dt 

and use conservation of the action 



tp ■ 



an 

dcost 



(87) 



S 



pdq 



over the period of motion for very slowly changing pa- 
rameters of the system [here h z (t)\. Since in the final 
state h z — > oo the spin precesses around the z axis with 
a constant value of s 2 (oo) = cos [0(oo)] that is related 
to the staying probability Pn by Eq. l)4*H*jl and the cor- 
responding action is simply 5(oo) = 27rcos#(oo), one 
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obtains 



P. 



N 



-s 



2tt 



oo, 



(89) 



where S is the action over the trajectory that starts from 
the metastable energy minimum that is on the verge of 
being unstable 



... = -(l-# 8 ) 1/a , 



= h 1/3 



0. (90) 



This trajectory can be found from the conservation of 
energy, 



- 2(1 - h 2 J z f' 2 t 



2K 



= l-3/^ 3 . (91) 



Analysis shows that in the range < h x < 3\/3/8 this 
trajectory encircles the z axis, thus the z axis can be 
used as the polar axis (s z = cos#, s x — sin # cost/?, s y = 
sin#sin<p) to compute S. In the overlapping range 1/8 < 
h x < 1 the trajectory encircles the x axis that can be 
used as the polar axis (s x — cos#, s z — — sin#cos</J, 
s y = sin 6 sin ip) . With these choices, cos 9 can be found 
numerically from Eq. (|91|l as a function of <p and the 
action S can be computed from Eq. I|88|l . The result for 
Pjv vs h x is shown in Fig. 1151 with asymptotes 



P 



N 



(3/2)hl f \ 



h x < 1 



0.544861 (1 - h x ) 3/2 , 1 - h x -C 1 



(92) 



VII. DISCUSSION 

The simplified model of interacting tunneling parti- 
cles that was considered above allows to quantify the in- 
fluence of interaction on the Landau-Zener-Stueckelberg 
staying probability P. It was done here by numeri- 
cally solving the problem for up to N = 200 interact- 
ing particles as well as using a number of analytical and 
half-analytical approaches, including the mean-field limit 
N -> oo. 

In accord with physical expectations and considera- 
tions of the exact levels of the system in Fig. [IJ,, the fer- 
romagnetic coupling tends to suppress LZS transitions 



to another bare energy level (i.e., the state on the other 
side of the energy barrier) , in agreement with Ref . 0- For 
N — > oo there is a critical value of the coupling above 
which the staying probability P does not go to zero in 
the slow-sweep limit e — > oo, in contrast with the stan- 
dard LZS case. For finite N the dependence of P on 
e is a sum of many exponentials with greatly differing 
relaxation rates that is very slow approaching zero (see 
Fig. I10[l . The same should be the case for more realistic 
interactions of the ferromagnetic type. 

The negative coupling in our model (that corresponds 
to the antiferromagnetic frustrating coupling) tends to 
boost the LZS transitions. In the limit N — > oo the 
staying probability P even turns to zero at finite values of 
the sweep rate while the effective sweep rate becomes an 
odd function of time with retardation in the resonance- 
crossing region (see Fig- EJ ■ O n the other hand, models 
with more realistic antiferromagnetic interactions such 
as the nearest-neighbor interaction have the energy-level 
scheme strongly differing from that shown in Fig. ^d. 
Preliminary results show that such interactions tend to 
hamper LZS transitions instead of boosting it, although 
not to such an extent as ferromagnetic interactions. 

Theoretically the model considered here is interesting 
since it maps on the problem of a single large spin S — 
N/2 and thus the mean- field limit N — > oo corresponds 
to the classical limit S — > oo for the large spin. It would 
be very interesting to study deviations from the mean- 
field solutions for large but finite N. These deviations 
can be very large, as can be seen in Fig. [HI in the region 
e ~ 0.11. How quantitatively well does the MFA work 
for model systems with more complicated interactions 
and for realistic systems remains unclear, and it is an 
interesting topic for further work. 
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